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Abstract —Beacon node placement, node-to-node measurement, 
and target node positioning are the three key steps for a 
localization process. However, compared with the other two steps, 
beacon node placement still lacks a comprehensive, systematic 
study in research literatures. To fill this gap, we address the 
Beacon Node Placment (BNP) problem that deploys beacon nodes 
for minimal localization error in this paper. BNP is difficult 
in that the localization error is determined by a complicated 
combination of factors, i.e., the localization error differing greatly 
under a different environment, with a different algorithm applied, 
or with a different type of beacon node used. In view of the 
hardness of BNP, we propose an approximate function to reduce 
time cost in localization error calculation, and also prove its 
time complexity and error bound. By approximation, a sub- 
optimal distribution of beacon nodes could be found within 
acceptable time cost for placement. In the experiment, we test our 
method and compare it with other node placement methods under 
various settings and environments. The experimental results show 
feasibility and effectiveness of our method in practice. 

I. Introduction 

Localization is a critical enabler for today’s context-aware 
applications, attracting tremendous research effort in recent 
years. Researchers have devised various approaches to improve 
localization accuracy, e.g., adopting more accurate signal mea¬ 
surement, using advanced techniques to alleviate measurement 
error, inventing new models to position target node, etc. For 
such approaches, beacon node placement is the prerequisite 
to performing localization. The localization error would differ 
with different distribution of beacon nodes. By careful place¬ 
ment of beacon nodes, localization accuracy can be improved, 
sometimes significantly. 

Despite its practical importance for improving localization 
accuracy, beacon node placement has not been thoroughly 
studied yet in research. Only some of existing works on node 
placement have considered the problem in the context of lo¬ 
calization in limited aspects. For example, the boundary effect 
on beacon node placement was discussed in 0; random, max 
and grid placement were compared in (4); optimal placement 
in camera network was studied in G). There still lacks a 
systematic study on beacon node placement for localization. 

In our experience of beacon node placement, we have 
found that the challenges lie in environmental context. The 
localization error is determined by complicated combination of 
factors. Modification on a single factor like adding or remov¬ 
ing beacon nodes, applying a different localization algorithm 
would vary localization error greatly. As a consequence, we 
usually cannot fully understand the cause and the effect of 
factor change on localization error. It makes establishing direct 


relationship from beacon node placement to localization error 
seems unapproachable. 

In this paper, as the first attempt, we set out to tackle the 
challenges, and try to find the optimal beacon node place¬ 
ment that has the minimal localization error with any given 
environmental context. In summary, we make the following 
contributions: 


• We model Beacon Node Placement (BNP) problem, and 
prove it is NP-hard. Therefore, we choose to find ap¬ 
proximate solution to achieve minimal localization error 
within acceptable time cost. 

• We propose an approximate approach that is orthogonal 
to environmental context, independent of measurement, 
localization algorithm, etc. We also provide the execution 
time complexity and localization error bound of the 
approach. 

• In experiment, we test and compare our method with 
other node placement methods under various settings 
and environments, such as 2910m 2 indoor floor and 
outdoor city-wide dataset. The experimental results show 
feasibility and effectiveness of our placement method. 


The rest of the paper is organized as follows. Section [II] 
formulates Beacon Node Placement problem, and prove its 
NP-hardness. Section [III] initially analyses the approximation 
on error calculation, and then propose several approximate 
techniques that are orthogonal to environmental context based 
on the analysis. Section IV discusses practical considerations 
in implementation, and synthesizes the approximate function 
by a combination of the approximate techniques. Section [V] 
evaluates our placement method and compares it with other 
methods under various settings. Section [VI| introduces the 
related work. Section IVIII summarizes our work. 


II. Problem Formulation 

The localization error is determined by a complicated 
combination of factors. It is unrealistic to give an explicit 
expression. Here, we use an abstract function to represent the 
calculation on localization error. 


e = f(M,A,I(B),l(B )) (1) 

where the localization error e is calculated by the input map 
M, the localization algorithm A, and the given beacon nodes 
B including their properties 1(B) and placement locations 
1(B). In this paper, we focus on the relationship between the 
localization error e and beacon nodes placement 1(B) in /(.), 



so as to find a distribution of beacon nodes with the minimal 
localization error. We formally define our target as follow. 

Problem 1 (Beacon Node Placement (BNP)). Given a map M, 
a localization algorithm A, and the information 1(B) about 
beacon nodes, the beacon nodes should be placed on locations 
1(B) that minimize the localization error 

Theorem 1 . BNP is NP-hard. 

Proof: (Sketch.) We give an instance of BNP as follows. 
Let M be a random, bounded area, and A be the trilateration 
algorithm. Assume that the set B of beacon nodes adopts the 
distance measurement model, the unit disk coverage model, 
and the non-sleep energy model. Then, M can be considered 
an infinite set IS of location points. Under the specified A 
and 1(B) , the localization error of every location point in IS 
can be computed independently to determine the localization 
error e. 

The above instance can be reduced to the set cover problem 
E2. The decision version of the set cover problem is NP- 
complete, and the optimization version of set cover problem 
is NP-hard 113 ]. In the decision version of set cover problem, 
every element in the universe should be checked if it is 
covered by selected subsets. By direct reduction of the error 
computing function of the above BNP instance to this function 
of checking coverage for every (sampled) location point in IS, 
we can prove that BNP is NP-hard. ■ 

III. Approximate /(.) 

A. Basic Idea 

In view of the hardness result for BNP, in this paper we try 
to employ an universal approximate approach that can deal 
with any map M, any localization algorithm A, and any given 
1(B) about beacon nodes. Our approach aims to synthesize a 
function /'(.) instead of /(.) to approximate the calculation 
of localization error e, so as to find an appropriate solution 
within acceptable time cost for BNP. 

Here, we first discuss the error of approximating /(.). We 
define the approximate error Ae = \e — e'\ = \f(x) — f'(x')\, 
where x = [M, 1(B), A, l(B)\. To further analyze the approx¬ 
imation error, we introduce Aft and Lf> to depict /'(.). Aft 
is used to describe the difference between /(.) and /'(.), and 
we have 

\f{x)-f(x)\<A r \x\ (2) 

where \x\ denotes a representative value of input. Lp is the 
Lipschitz constant to describe the smoothness of /(.) between 
any two different inputs x and x', and we have 

\f(x) - f(x')\ < Lf\x - x'\ (3) 

where \x — x'\ represents the difference between x and x'. 
Here, we use an abstract representation of \x\ and \x — x'\ 
since their values may be influenced by many hidden factors 
that are hard to infer in an universal approach. And to describe 
the impact of different factors, we introduce weight vector 
W G R m , having 

e = f(x) = f{{xj x Wj | j = 1,2,to}) (4) 

'The localization error e could refer to arithmetic average error, geometric 
average error, median error from a set of results, etc. The techniques discussed 
later can be applied no matter which specific meaning it has. 


TABLE I 

Summary of notations 


Notation 

Description 

M 

the input map for beacon node placement 

M s 

the signal map generated by beacon nodes 

A 

the localization algorithm applied 

B 

the set of beacon nodes 

~XB) 

the information about beacon nodes including mea¬ 
surement model, energy model, coverage model, etc 

“P) 

the output location for beacon node placement 

/(■) 

denoting the function /(M, I(B ), A, 1(B)) for short 

e = /(•) 

the localization error e for all locations in M can be 
calculated by the function /(.) 

e' = /'(•) 

the approximate error e' is calculated by the approx¬ 
imate function /'(.) 

e opt 

the minimal localization error can be achieved in the 
selected area 

\> 

1 

II 

< 

the absolute difference between e and e' 

A f , 

the difference factor between /(.) and /'(.) 

L f , 

the Lipschitz constant for the function /(.) 

w 

the weight vector for the inputs of /(.) 

V 

a location point in map M 

<?(•) 

denoting the relationship between localization error 
and distribution of beacon nodes 

~hU) 

denoting the signal distribution around a beacon node 

Coll 

the collection of signal points related to a beacon node 

Pace 

user-specified acceptable calculation time 

^sBacc 

user-specified acceptable localization error 

IsPacc 

the sampling interval calculated by g(AE acc ) 


The function /'(.) could be composed by a sequence of 
approximate techniques. We use < 

to denote these techniques, and let e ' denote the error change 
when applying f'(x) following < f[(x), / 2 O),/•_ ± (x) >. 
We have 

e i = fi({xj x w i,j I * = 1) 2, ••., n, j = 1,2,..., to}) (5) 
Combine Eq. (j2ji and Eq. <[3j>, we get 

Ae =|e-£”=U'l 

= \f(x)-E7=ifl(x')\ 

<\f(x) - Ya=i nix) 1 +£” =1 mix) ~ nix’) \ 

< xr=i ET=ii A fi\ x i w i,i\+ L fi\^iWij\) 

( 6 ) 

As can be seen from Eq. @ 1 , in general, the approximate 
error Ae is bounded by the sum of impact of each approximate 
technique and each input factor. For the two powers A ft and 
Lp, if either has a large value, the impact of input factors 
would be amplified and thus Ae would become large. So 
the design of approximation should be carefully made to 
keep the approximate function /'(.) smooth and very close 
to /(.). However, we should mention the smoothness of /(.) 
self here. If the applied localization algorithm cannot work 
stably under various environments (or various input factors), 
then /(.) self would be unstable, which will make Lp and 
Ae have a large value. But we believe it is less likely 
to happen when applying those frequently-used localization 
algorithms, otherwise they won’t be adopted in practical use. 
As for the impact of input factors, stated plainly, the one 
Xi with bigger weight Wij would have greater influence. 
Nevertheless, for our approximate approach, it is impossible 
to figure out the exact input factors and their weights which 


























Fig. 1. A Sampling example. In this example, location point pi, p2, ..., 
pig are uniformly sampled in the hexagon; 3 beacon nodes are planed to be 
placed in this hexagon. By sampling, there are altogether Cf Q combinations 
for beacon node placement; for each combination, the worker measures the 
signal in these 19 points. In the figure, 3 beacon node are placed on p$, pn 
and pi4 respectively, and the worker is at pg. 


depend on the specific localization scene. For this reason, we 
shall consider every location point in M as equally important 
in approximation so as to avoid deviation on input factors. 
Later we will show how to design techniques on approximation 
according to this discussion. 


B. Techniques on Approximation 

We propose several techniques on approximating /(.) in 
this section. These include Sampling, Memorization, Skipping 
and Interpolation. All of these techniques are proposed for 
reducing the execution time cost for BNP, and can be applied 
individually. 

1) Sampling: For some localization scenes, it is impossible 
to calculate out the error on every location point in M within 
an accepted time cost. A natural way is to efficiently sample 
useful location points in M to approximate the localization 
error e. As discussed, we are not able to infer the exact input 
factors and their weights in the universal approach. So we 
apply uniform sampling on M to treat every location point 
as equally important. The approximate function /'(.) with the 
sampling technique applied would have a reduced calculation 
time compared to directly calculate e by the function /(.). 
We illustrate a sampling example in Fig. [l] In this example, 
a hexagon on the map is sampled with uniformly distributed 
location points. 

2) Memorization: Usually, the coverage range of a beacon 
node is limited. Therefore, we can memorize the calculation 
results of a selected area to infer the results of other areas. 
More specifically, when calculating the localization error with 
a given distribution of beacon nodes, we can look up the mem¬ 
orized results for the same or similar beacon node distribution 
as an approximation to reduce the calculation time. As shown 
in Fig. [2j 3 beacon nodes are located at point p\, p2, and p% 
separately inside a hexagon. If the memorized results of the 
hexagon (with dashed line) have 3 beacon nodes with their 
distribution \p[ -pi\ + |p' -p 2 \ + \p' k - P:i | < T mem , where 
Tmem is a threshold distance value set up, then the error of 
distribution p\, p 2 and ps can be approximated by the result 
on the distribution pi , pj and p k . 

3) Skipping: A bad distribution of beacon nodes would not 
be a good approximation of /(.), and thus can be ignored in 
error calculation to reduce the time cost. Intuitively, we can 
infer a distribution to be a bad o ne based on previous results. 
As discussed in Section |III- A we assume that frequently- 
used localization algorithms are stable ones, otherwise they 
would not be used in a practical manner. Therefore, for a 



Fig. 2. A Memorization example. Suppose that the calculation results of 
beacon node placement in the hexagon with dashed line have already been 
memorized. When deploying beacon nodes in the hexagon with solid line, 
similar node distributions can be looked up for in the memorized results. In 
the figure, we translate pi, pj , p k to p', p '-, p' k by the vector pfpi respectively. 
If I p\-p 1 1 +1 p’. -p 2 1 +1 p' k -P3 1 < Tmem, we consider that the distribution 
of pi, p2 and p3 could be approximated by pi, pj and p k . 



Fig. 3. A Skipping example. Suppose that some distributions of beacon 
nodes are shown as bad ones during previous calculation. Then, when beacon 
nodes are deployed with the distribution similar to bad ones, we can infer 
that these nodes are also in a bad arrangement. In the figure, the set of points 
nearest in different directions are found for pi , p2 and p3 respectively. For 
any combination of pi G {P4,P5,P6,P7}, Pj G {P8,P9,Pio,Pn} and 
Pk G {P12,P13,P14,P15}, if I Pi ~Pi\ + | Pj ~P2\ + | Pk P3 1 < T skp 
and the combination of pi , pj and p k has unacceptable localization error, we 
consider that the calculation on pi, p2, and p3 can be skipped. 


specific distribution of beacon nodes, if all its similar distri¬ 
butions have unbefitting error, we believe that this distribution 
could be skipped without calculation. We also set a distance 
threshold value T s p p for Skipping. Generally speaking, we 
have T s p p > T mem since the result for reuse should be more 
accurate than the result for inferring its bad or good. Take 
Fig. [3] for example. For any distribution pi G {P 4 ,P 5 ,P 6 ,P 7 }, 
Pj G {P8,P9,Pw,Pn}, Pk G {Pi2,Pi3,Pi4,Pi5}, if we have 
|Pi ~Pi\ + |Pj -P 2 I + \Pk ~P 3 \ < T skp and the distribution 
of pi , pj, and pk is considered as an unaccepted one for 
approximation, we choose to skip the error calculation on the 
distribution of pi, P 2 , and £> 3 . 

4) Interpolation: The function /(.) self could be involved 
with complex error calculation. Either /(.) has a non-linear, 
complex form which may take a lot of time on calculation, 
or it cannot be directly derived in theory which needs a 
numerical simulation approach. To deal with the case, we 
use polynomials to approximate /(.) based on its form or 
its numerical simulation results. For example, e x can be 
approximated by its Taylor series Yli= o x V^-’ if 3c having 
I YTi =0 x l /i\ — e x \ < Ae where Ae is a defined limited error 
on approximation. 











Algorithm 1: ErrorOnDistribution 
Input: Input map M, the localization algorithm A, the 
set of beacon nodes B, information 1(B) about 
beacon nodes 

Output: The relationship g(.) between localization error 
and distribution of beacon nodes 

1 for j = 1, 2 , ...n do 

2 Randomly assign locations p^ i = 1, 2, \B\, to 
beacon nodes B; 

3 Randomly generate an offset distance A pp, 

4 Randomly move beacon node with 

Eifi \p'i -Pi\/\ B \ = A ft; 

5 Calculate the difference error Ae between 
{Pi,i = 1) 2,|.B|} and {pi = 1,2, \B\}; 

6 Store (Api, Ae^); 

7 end 

8 Fit Ap = g(Ae) by {(Ap*, Ae*)|i = 1,2, 

9 return g(.); 


Algorithm 2: SelectedArea 

Input: Acceptable calculation time T acc and localization 
error AE acc , the relationship p(.) between 
localization error and distribution of beacon 
nodes, input map M, the localization algorithm 
A, the set of beacon nodes B , information 1(B) 
about beacon nodes 
Output: The selected area Area 

1 Calculate the offset distance A P acc = g(AE acc ); 

2 Set the sampling density d s = l/AP a 2 cc point/size; 

3 Set the beacon node density d 5 = |T?|/|M| node j size; 

4 Time the execution of /'(.) as t; 

5 Solve the area size S by max(S) s.t. C^ b • t < T acc ; 

6 Select an area Area (such as a regular hexagon) of size 
S; 

7 return Area; 


IV. The Proposed BNP Algorithm 
A. Practical Consideration 

With the NP-hardness of BNP, it is infeasible to directly 
calculate on /(.) for optimal beacon node placement. In¬ 
stead, we apply the techniques discussed in Section |III-B| to 
approximate /(.). Nevertheless, to combine these techniques 
in the system, three critical issues should be considered for 
us in implementation. The first is to let the approximate 
function /'(.) execute within acceptable time. The second is 
to characterize localization environment on the execution of 
/'(.). The third is to consider the boundary of map for BNP. 
Next we address these three issues respectively. 

1) Strategy on Error-Time Trade-off: As discussed in Sec¬ 
tion |III-B2| the distributed pattern of beacon nodes in a se¬ 


lected area can be memorized to infer beacon node placement 
in other areas. At the same time, when applying memorization, 
a trade-off in area size inevitably occurs: either to explore a 
larger area for memorization to provide more accurate approx¬ 
imation of /(.), or a smaller area to reduce the calculation 
time of /'(.). To deal with this trade-off, we should select an 


Algorithm 3: ModelingNode 
Input: Input map M, the set of beacon nodes B , 
information 1(B) about beacon nodes 
Output: the collection Coll of signal points, the fitting 
expression h(.) of signal related to a beacon 
node 

1 User places a beacon node b on a point p in M; 

2 Gather the signal points around b to the collection Coll; 

3 Fit the relation of signal points to b as the expression 
h{.)- 

4 return Coll and h(.); 


Algorithm 4: Modeling Area 

Input: Input map M, the set of beacon nodes B , 
information 1(B) about beacon nodes, the 
collection Coll of signal points, the fitting 
expression h(.) of signal related to a beacon node 
Output: The signal map M s 

1 for Each signal point p in M do 

2 if Similar points found in Coll then 

3 Generate signal by the expected value of similar 
points; 

4 else 

5 | Generate signal by h(.); 

6 end 

7 Record the generated signal to M s ; 

8 end 

9 return M s ; 


area of proper size to balance between localization error and 
calculation time. 

To select a proper sized area, we first try to determine 
the relationship between localization error and distribution of 
beacon nodes by random assignment. As illustrated in Fig. 
4(b)| three beacon nodes are randomly assigned with location 
Pi , pj and pk . We set an offset distance Ap, and randomly 
move the beacon nodes to location p', p’- and p' k having 

J2q=i I Pq — P q \/3 = Ap. Then, for these two beacon node 
distributions, we calculate their difference value, Ae, on error. 
By sampling a group of such beacon node distributions, we 
have {(A p q , Ae q )\q = 1, 2, ...,n}. We use polynomial fitting 
to approximate the relationship between Ap and Ae, formally 
as Ap = g( Ae). We summarize the process in Algorithm [T| 
— ErrorOnDistribution. 

Based on g(.). Algorithm^- SelectedArea - gives a sketch 
on determining the size of selected area for memorization. 
SelectedArea takes the acceptable calculation time T acc and 
acceptable localization error A E acc as the user-specified input. 
Within the limit of T acc , the approximate function /'(.) has 
|/'(.) — /(.)| < A E acc . In process, SelectedArea uses the 
sampling density d s and the beacon node density d 5 to estimate 
the area size S by max(S), s.t. • t < T acc , that is, to 

select an area of proper size by considering both T acc and 
A E acc with respect to the execution time of /'(.). 

In addition, we use a • C^ b • t < A E acc with an 
approximate ratio 0 < a < 1 in practical calculation, since 
























techniques on approximation are ap plied, i.e., some execution 
of /'(.) are skipped (Section III-B3). We omit the detail about 
a due to it is not a primary concern in our paper. 

2) Strategy on Characterization of Localization Environ¬ 
ment: Naturally, we can place beacon nodes in different 
distributions to characterize the environment, and thus to find 
a relatively lower localization error e!. However, it would 
cost abundant labor in placement with repeatedly adjusting 
the position of beacon nodes and calculating their localization 
error distributions. Also, freely moving beacon nodes are 
impossible in some localization scenes, i.e., beacon nodes need 
electric plug. So we believe it is not an applicable way for our 
general approach of BNP. 

Instead, we apply a modeling approach based on mea¬ 
surement results. First, we gather the strength (or sometimes 
the direction) of signal points scattered over an area or a 
straight line from a beacon node as a collection Coll , and 
apply Interpolation (Section III-B4) to fit the expression h(.) 
of signal strength (or direction) related to a beacon node. 
We summarize these course of actions in Algorithm [3] - 
ModelingNode. Then, when modeling the whole map, for 
every point, we search in the collection Coll for similar ones, 
and generate signal by the expected value of these similar 
ones; otherwise, if no similar point found, we generate a signal 
by the fitting expression h(.). Algorithm]^ - ModelingArea - 
depicts the process of modeling beacon node placement. Later 
in the experiment (Section [V]), we also describe the process 
of ModelingArea in practical environment. 

3) Strategy on The Boundary of Map: In practical deploy¬ 
ment, the input map M should be considered with boundary. 
When there is an existing boundary, we need to address 
two cases in BNP separately. In one, beacon nodes may be 
placed outside of the boundary. Therein, BNP can be regarded 
as working on the unbounded area, and we only need to 
place beacon nodes to cover M. The other is that beacon 
nodes cannot be placed out of the boundary, i.e., beacon 
nodes cannot be placed outside in the second floor of a 
building. In this case, we search in the memorized results 
for the beacon node distribution with two criteria: all beacon 
nodes are located inside the boundary; the localization error 
should be minimized on this occasion. Otherwise if no such 
distribution is found, we search for the best memorized result 
with minimal localization error, and move each outer beacon 
node of this result inside following the shortest distance. This 
process on the boundary of map is also described in Algorithm 
[5] — ProcessingBoundary. 


B. The Algorithm and Complexity 

Combining all the discussions above, we can synthesize 
the approximate function /'(.) and find a sub-optimal beacon 
node placement by Algorithm [ 6 ] The main idea behind it is to 
divide-and-conquer the calculation of /(.) on M, by which we 
can lower down the calculation time to a user-specified T acc \ 
place beacon nodes to locations 1(B) with the localization 
error kept at most A E acc (or yAE acc addressed later) greater 
than the error by the optimal beacon node placement. 

To clearly understand the process of synthesization in Al¬ 
gorithm [ 6 ] we explain it with demonstrations in Fig. [4] In 
line 1 , we collect signals from a beacon node at scattered 
location points as collection coll and fit the expression h(.) 


Algorithm 5: ProcessingBoundary 
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Input: Input map M, the set of beacon nodes B , 

information 1(B) about beacon nodes, memorized 
result memo , the output location 1(B) 

Output: The output location 1(B) 

Let db est be the distribution with minimal localization 
error in memo ; 

if Beacon nodes are allowed to be placed outside M 

then 

| Write 1(B) according to db est ; 

else 

Look up memo for the distribution d with beacon 
nodes inside the boundary and minimal localization 
error; 

if d 7 ^ Null then 
| Write 1(B) according to d\ 

else 

Write 1(B) by moving each outer beacon node of 
dbest inside following the shortest distance; 

end 

end 

return 1(B ); 


by Algorithm ModelingNode (Fig. |4(a)| ). In line 2, we model 
the relationship g(.) between the offset distance A p and the 
difference error Ae by Algorithm ErrorOnDistribution (Fig. 
4(b)| ). In line 3, taking user-specified acceptable calculation 
time T acc and localization error A E acc as input, we select 
a circle area with its size con strain ed by T acc and A E acc 
in Algorithm SelectedArea (Fig. 4(c)| ). In line 4-32, we apply 
sampling on area. At the beginning of sampling process, we 
firstly compute the acceptable sampling interval A P a cc , and 
select a initial interval inti (inti > A P acc ) (Fig. 4(d) ). Then, 
during sampling process, we examine each combination of 
beacon node distribution (line 12-26). For a distribution, if all 
the similar distributions are considered to be bad ones, we skip 
its calculation (Fig. 4(e)| ); else we calculate the locali zation 
error on this distribution and memorize the result (Fig. |4(f)| ). 
Besides, in line 27-39, we use polynomials to approximate 
the error calculation after obtaining a set of results. Next, in 
line 33-41, we apply the best distribution found in area as a 
pattern to the rest of map M — area. In it, we consider two 
problems in practical deployment. The first is to deal with 
the boundary of map M. When a beacon node is not allowed 
to be placed out of the bou ndary, we either move this node 
inside or just remove it (Fig. 4(g)| ). The other is to adjust the 
distribution pattern of beacon nodes and the area based on 
practical environment (Fig. |4(h)| ). 


In brief, the synthesization process looks up for the best 
distribution in a selected area, and applies it as a pattern 
to the rest of the map. A hypothesis behind it is that the 
error distribution of the selected area is the same as other 
areas in the map. Thus, we define an ideal case: the error 
distribution in each circle (or regular hexagon) is the same 
(namely the function g(.) can be applied globally); beacon 
nodes can be placed out of the boundary. Let e opt denote the 
minimal localization error can be achieved in the selected area. 
We have Theorem [2] on time complexity and localization error 




























(a) A beacon node is placed on the point p± 3 . (b) At first, three beacon nodes are randomly assigned with (c) Given a user-specified acceptable calcula- 
The worker gathers signals from the beacon location p^, pj and p Then, we set an offset distance A p, tion time T acc and localization error A E acc , 
node by sampling at location point pi, p 2 , and randomly move these three beacon nodes to location we compute the offset distance A P ac c = 
..., P 25 respectively. Then, we get the col- p'., p'. and p' k having Ylq=i I Pq ~ Pq I/ 3 = Ap. Next, g(AE acc ) which means the sampling den- 
lection Coll, and can fit the expression h(.) f 0 r these two beacon node distribution, we calculate their sity should be at least as low as 1/AP^ CC 
(Algorithm ModelingNode). For reducing the difference error value Ae. By sampling a group of such point /size', then, as shown in the figure, 
labour cost, this result can be directly used to beacon node distributions {(Ap g , Ae q )\q = 1, 2, n}, we generate a circle area with its size 
represent the localization environment such we apply polynomial fitting to approximate the relationship |^rea| constrained by T acc . (Algorithm Se¬ 
as modeling (A p, Ae) in Fig. |4(b)| Ap = g(Ae). (Algorithm ErrorOnDistribution) lectedArea) 



(d) We sample on area with a initial interval (e) We apply finer granularity sampling with interval (f) We choose to calculate the localization error 
inti (inti > AP acc ). In the figure, it has inti /2 (inti /2 > AP acc ). For each distribution of beacon e for a distribution of beacon nodes if it cannot 
total 13 sample points. If 3 beacon nodes nodes, we look up for similar ones in the memorized be skipped. Before calculation, we take the dis- 
are planned to be placed in area, there results. And if there exist similar ones with unacceptable tribution of beacon nodes in area as a pattern 
are altogether Cf 3 combinations for beacon localization error, we believe current combination could be and apply it to neighbor areas. Then, based on 
node placement, and for each combination skipped without calculation. In the figure, we search the the signals gathered from beacon nodes in area 
the worker measures the signals from these set of points nearest in different directions for pi, p 2 and and neighbor areas at each sample point, we 
3 beacon nodes at each sample point. Then, p 3 respectively. If all the similar distributions out of these calculate the localization error e and memorize 
we calculate the localization error e for each sets are considered as bad ones, we skip the calculation of the result, 
combination and memorize the result. localization error for pi, p 2 and p 3 (also shown in Fig.[3j. 




(g) After the best distribution of beacon nodes with minimal localization error is found (h) Unlike the ideal case in Fig. |4(g)| in practical environment, due 
in | area\, we apply it as a pattern to the rest of the map (namely M — area). When to the reason of interference, block; etc., signal from beacon nodes 
an area encounters the boundary and beacon nodes are not allowed to be placed out may differ significantly. Therefore, we adjust the distribution pattern 
of the boundary, we compute the overlapped size of this area with M, and look up in of beacon nodes and the area based on circumstances. In the figure, 
the memorized results for applicable distribution with minimal localization error on this an obstacle blocks Line Of Sight (LOS) transmission of signals from 
occasion. Otherwise if no such distribution found, we simply move outer beacon nodes beacon nodes. Considering the area with three beacon nodes located 
of the applied pattern inside following the shortest distance. In the figure, beacon nodes at p±, p2, P 3 , suppose that a beacon node is placed at its central 
at P 3 , P 4 , P 5 and pq are removed since the area they located have small overlapped size point, we use the contour line of — 70dBm signal from this node 
with the map M, i.e., the circle with p 5 and p 6 has the overlapped size less than | of its as the termination condition to limit the coverage range of the area, 
whole size; the beacon node at pi is moved to p [ with applicable distribution is found Then, we adjust these three beacon nodes located at p±, p2, P3 to 
in the memorized result; the beacon node at p2 is moved inside to p ' 2 with the shortest the location of p' x , p' 2 , p ' 3 by bilinear interpolation with respect to 
distance since no applicable distribution is found. (Algorithm ProcessingBoundary) the change of the area. 


Fig. 4. Synthesize /'(.) 



















for the synthesization process. 


Theorem 2. In the ideal case , the synthesization process 
(Algorithm [b]) executes with time complexity of 0(T acc ) and 
generates the distribution of beacon nodes with localization 
error e - e opt < yAE acc , 7 > 0. 


Proof: In the selected area, we have the total calculation 
time • t < T acc (Line 5, Algorithm SelectedArea). The 

best distribution found in the selected area can be directly 
applied to other areas. Therefore, the synthesization process 
has time complexity 0(T acc ). 

Let P opt be the best distribution of beacon nodes found, 
having the minimal localization error e opt . With the sampling 
interval set to A P acc (Line 1, Algorithm SelectedArea ), we 
can find a distribution P with localization error e, having 


IP - Pr 


opt | _ 


< 


\/2 p 

2 1 acc’ 


It has the error e > 


_ &opt 


if the 


approximate function /'(.) does not change the monotonicity 
of /(.). As can be seen, the approximate techniques (Sampling, 
Memorization, Skipping, and Interpolation) applied in the syn¬ 
thesization process do not change the original monotonicity of 
error distribution. (For Interpolation, if it has low error.) Thus, 
we can infer that the approximate function /'(.) is monotone 
increasing around P op t- (Namely, |P— P op t\ °c e—e op t •) Recall 
that A P acc = g(AE acc ). For |P - P opt \ < ^P acc , we have 
e &opt E 7AP acc , 7 — 0. ® 

In a practical environment, it is complicated to determine 
7 in Theorem [2] since the value of 7 depends on the given 
map, localization algorithm, and the information about beacon 
nodes. As discussed in Section III the frequently-used local¬ 
ization algorithms are less likely to be unstable in smoothness, 
otherwise they would not be used in practical. Thus, we believe 
that drastic variations on error distribution are unlikely to 
happen, and 7 usually has a small value. 


V. Experimental Analysis 

In this section, we evaluate our beacon node placement 
method in actual environments. First, we change the size of the 
selected area and the sampling interval to assess our method 
on execution time and localization error. Then, we compare 
our method with several other placement methods in indoor 
environment. Finally, we also experiment on a large scale, 
outdoor real-world dataset. 


A. Experiment Setting 

The experiment setting is as follows. 

Measurement and Localization Algorithm: There exist nu¬ 
merous studies on measurement and algorithm for localization. 
Among these studies, WiFi based localization has attracted 
tremendous attention for its wide availability and no extra 
deployment cost in recent years. Most of existing WiFi based 
localization algorithms can be divided into two categories: 
Model-based and Fingerprinting-based. In our experiment, we 
consider WiFi measurement and select EZPerfect 003 and 
RADAR 0, which are representative of model-based and 
fingerprinting-based algorithm respectively. 

• EZPerfect trains the parameters of the log-distance path 
loss model by sampling signals at selected locations, 
and apply Trilateration or Multilateration to estimate the 
location of target points. 
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Algorithm 6: Beacon Node Placement 
Input: Input map M, the localization algorithm A, the 
set of beacon nodes B , information 1(B) about 
beacon nodes, user-specified acceptable 
calculation time T acc and localization error A E acc 
Output: The output location 1(B) of beacon nodes 
[coll, h(.)} = Modeling Node(M, B,I(B))\ 
g(.) = ErrorOnDistribution(M, A, B,I(B )); 
area = 

S elected Area(M, A, B, 1(B), g(.),T acc , AE acc )\ 
n = |area| • |L?|/|M|; //Number of nodes to be placed 
A Pace — g(A E acc )\ I/Acceptable sampling interval 
User set the interval inti ( inti > AP acc ); 

Set the error e = oc; 

Set the location 1(B) = 0; 
while inti > AP acc do 
//Apply Sampling 

m = \area\/intl 2 \ //Number of points to be sampled 
foreach Beacon node distribution {pi\i = 1,2, ...,n} 
out of all C^ combinations do 

Modeling Area(M, B, 1(B),Coll, h(.)); 

//Apply Skipping 

if {pi\i = 1,2,... , n} can be skipped then 
| Continue; 
else 

Calculate the localization error e! by 
sampling (namely using /'(.)); 

if e! < e then 

e = e!\ 

Record the locations of beacon nodes to 

KB ); 

end 

//Apply Memorization 

Memorize {pfi = 1,2, ...,n} and e! to 

memo ; 

end 

end 

HApply Interpolation 

if the error calculation can be interpolated then 
Use polynomials to approximate the error 
calculation (namely approximate /'(.)); 

end 

inti = inti/ 2; 

end 

Split the remaining area M — area to 
{areai \ i = 1,2, 

Let db es t be the distribution with minimal localization 
error in memo ; 

foreach areai in {areai | z = 1, 2, ...,q} do 
if areai contains boundary then 
1(B) = 

ProcessingBoundary(M, B, I(B), l (B), memo ); 

else 

| Write 1(B) according to 

end 

end 

return 1(B ); 













96 m 



Fig. 5. At the beginning of the experiments, we randomly deploy some China Mobile M601 phones, with WiFi hotspot turned on, as beacon nodes. Then 
we use an AmigoBot with a phone placed on it traveling around the floor to gather signals from these beacon nodes. 


TABLE II 

The execution time of our placement method (Algorithm[6) on 
EZPefect in indoor environment 


EZPerfect 

lm 

2m 

3m 

4m 

5m 

S 

5.36 x 
10 51 (e.) 

5.27 x 
10 39 (e.) 

4.34 X 
10 32 (e.) 

3.57 X 
10 27 (e.) 

2.84 X 
10 24 (e.) 

is 

1.62 x 
10 n (e.) 

1.52 x 
10 8 (e.) 

2.52 x 
10 6 (e.) 

2.21 x 
10 5 (e.) 

1.78 x 
10 4 (e.) 

1S 

1.01 x 
10 5 (e.) 

1631.23 

125.96 

25.18 

12.08 

1 S' 

16 D 

1.09 x 
10 4 (e.) 

257.35 

15.07 

8.15 

0.75 

25 ^ 

3.02 x 
10 3 (e.) 

45.14 

6.14 

0.56 

0.25 


1 All values are in seconds, or s. The estimated values are with the 
suffix ’(e.)’. The selected area size and sampling interval vary from 
S to S, and Ira to 5ra respectively. 


• RADAR collects fingerprints of signals at known locations 
to establish a fingerprint database, and then determines 
the location of a target point by averaging the locations 
of these nearest fingerprints found in the database. 

Beacon Node Placement Method: For comparison, we im¬ 
plement four other placement methods, including two random 
ones and two deterministic ones. 

• Random optionally selects location points for beacon 
node placement. In the experiment, we generate distri¬ 
bution of beacon nodes for 3 times, and show the best 
result among them. 

. RKC GD assigns beacon nodes to location points that 
are randomly selected from the near-optimal hitting sets 
for k-covering area. We also generate beacon node distri¬ 
bution for 3 times and show the best result among them. 

• Uniform places beacon nodes at regular intervals. 

• CERACC |T] deterministically assigns beacon nodes to 
the lenses of slices in triangle lattice pattern for k- 
covering area. 

Localization Error: We use the following four error repre¬ 
sentations as the metric to evaluate the quality of beacon node 
placement method. 

• Arithmetic mean error (ari.): e ar i m = ^ Yl7=i G 

• Geometric mean error (geo.): e geo . = VTlILi G 

• Median error (med.): Let e\ < c^ < ... < e n . If n %2 == 
0, then e med . = en±i; otherwise e med . = (e- +en +1 )/2. 


TABLE III 

The execution time of our placement method (Algorithm|6J on 

RADAR IN INDOOR ENVIRONMENT 


RADAR 

lm 

2m 

3m 

4m 

5m 

S 

4.65 x 
10 51 (e.) 

4.30 x 
10 39 (e.) 

3.49 X 
10 32 (e.) 

8.20 X 
10 27 (e.) 

2.75 X 
10 24 (e.) 

i s 

1.15 x 
10 n (e.) 

1.16 x 
10 8 (e.) 

1.79 x 
10 6 (e.) 

1.73 x 
10 5 (e.) 

1.33 x 
10 4 (e.) 

is 

7.45 x 
10 4 (e.) 

1229.61 

97.44 

21.77 

9.38 

1 S' 

16 D 

8.37 x 
10 3 (e.) 

203.07 

13.08 

6.01 

0.54 

25 ^ 

2.35 x 
10 3 (e.) 

34.94 

4.92 

0.41 

0.19 


• Proportion of abnormal errors (abn.): Let bi =» 1 if e* > 
2e ari ', else set bi = 0. Then p abn . = bi x 100%. 

Experimental Environment: We experiment in two localiza¬ 
tion environments, with their difference on the type of beacon 
node, signal transmission, and scale. 

• Indoor Environment : We conduct the indoor experiments 
in the floor, shown in Fig. [5j with its area size S = 
2910m 2 . In experiments, we use a total of 20 mobile 
phones as beacon nodes to create WiFi hotspots. At the 
beginning of the experiments, we control an AmigoBot 
with a mobile phone walking around the floor, gathering 
signals at location points. Then, we divide the gathered 
signals into several collections colli , C 0 II 2 , coll n by 
distinguishing number of walls from the WiFi hotspot to 
the location point receiving signal, and fit the expression 

^ 2 (•)’ h n (.) (Algorithm ModelingNode). 

• Outdoor Environment'. We directly use the MetroFi 
dataset l22l . It involves 72 access points and samples 
signals at over 200,000 location points in a city-wide 
area. These location points are taken as candidate loca¬ 
tions for beacon node placement, and the access points 
are considered as the target nodes for localization in 
experiments. We take the whole dataset as collection coll , 
and fit the expression h(.). 

All the placement methods and localization algorithms were 
implemented with VC++. All the calculation were running on 
a Windows 7 machine with 2.3GHZ Intel Core i7 CPU and 
8GB RAM. 












































































TABLE IV 

The localization error ol our placement method (Algorithm^ on EZPefect in indoor environment 


EZPerfect 

2m 

3m 

4m 

5m 

an. 

geo. 

med. 

abn. 

an. 

geo. 

med. 

abn. 

an. 

geo. 

med. 

abn. 

an. 

geo. 

med. 

abn. 

I s 

2.85 

2.28 

2.47 

6.82% 

3.05 

2.45 

2.59 

6.58% 

3.12 

2.48 

2.68 

7.06% 

3.95 

3.15 

2.97 

7.40% 

I Q 

16 D 

2.94 

2.33 

2.57 

7.08% 

3.88 

3.01 

3.03 

8.82% 

4.05 

3.24 

3.09 

9.09% 

4.21 

3.31 

3.26 

10.049o 

1 o 

_25L_ 

3.10 

2.51 

2.68 

8.81% 

4.11 

3.41 

3.08 

9.28% 

4.29 

3.49 

3.80 

11.789 

0 4.50 

3.66 

3.89 

12.529o 


1 All values of error are in meters or m, except abn. using percentage. 


TABLE V 

The localization error of our placement method (Algorithm[6J on RADAR IN indoor environment 


RADAR 

2m 

3m 

4m 

5m 

an. 

geo. 

med. 

abn. 

an. 

geo. 

med. 

abn. 

an. 

geo. 

med. 

abn. 

an. 

geo. 

med. 

abn. 

A 

3.09 

2.31 

3.01 

0% 

3.18 

2.58 

2.82 

0% 

3.23 

2.72 

2.90 

0% 

3.50 

3.06 

2.77 

4.92% 

16 ^ 

3.13 

2.70 

3.27 

3.70% 

3.49 

2.96 

2.94 

1.72% 

3.51 

3.05 

3.45 

4.01% 

3.55 

3.14 

3.18 

5.62% 

— S' 

_2522_ 

3.19 

2.76 

3.05 

0% 

3.53 

3.03 

3.29 

4.12% 

3.82 

3.12 

3.56 

4.76% 

4.09 

3.25 

3.79 

7.41% 



(a) Random (b) RKC 
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TABLE VI 

The localization error of the referenced methods on 
EZPerfect in indoor environment 


EZPerfect 

ari. 

geo. 

med. 

abn. 

Random 

6.26 

5.16 

5.25 

15.09% 

RKC 

4.42 

3.56 

3.77 

8.08% 

Uniform 

4.66 

3.77 

3.86 

6.56% 

CERACC 

5.91 

4.63 

4.12 

8.06% 


TABLE VII 

The localization error of the referenced methods on RADAR 

IN INDOOR ENVIRONMENT 


(c) Uniform 


(d) CERACC 


Fig. 6. Beacon Node Placement 


RADAR 

ari. 

geo. 

med. 

abn. 

Random 

5.52 

4.41 

4.89 

12.32% 

RKC 

4.14 

3.29 

3.74 

9.57% 

Uniform 

4.11 

3.28 

3.81 

7.60% 

CERACC 

4.95 

4.21 

4.01 

8.58% 


B. The Impact of Selected Area Size and Sampling Interval 
on Execution time 

Here, we introduce the experimental results of our beacon 
nodes placement methods on execution time in an indoor 
environment. In experiments, we vary the size of selected 
area from S to — S', and sampling interval from 1 m to 5m, 
timing the execution of our placement method on localization 
algorithm EZPerfect and RADAR. The results of execution 
time on EZPerfect and RADAR are shown in Table [II] and [III] 
respectively (with some non-computable items estimated by 
Cg^ b • t in Algorithm SelectedArea). As can be seen from 
the results, by applying techniques (Sampling, Memorization, 
Skipping, and Intepolation) on approximation, we can largely 
reduce the execution time on finding beacon node distribution 
for placement, i.e., the execution time of our placement 
method with selected area size |S and sampling interval 2m 
(15, 2m) on EZPerfect reduced by a factor of 3.29 x 10 48 
compared to the case of (5,1m). This experimental result 
shows the effectiveness of the approximate techniques applied. 
Besides, we do not verify the effectiveness of each technique 
independently here since it is not the major concern in our 
evaluation. 

C. The Impact of Selected Area Size and Sampling Interval 
on Localization Error 

Next, we focus on the impact of varying the size of a 
selected area and sampling interval on localization error. In 


the experiment, we target finding the beacon node distribution 
with minimal arithmetic mean error. (Besides, we can get 
similar results using other error representations. Due to space 
limits, we omit them here.) The results on EZPerfect and 
RADAR are shown in Table IV and [V] respectively. As can 
be seen from both of these two table, 1) the localization 
error increases as the selected area size changes from ^5 
to ^5; 2) the localization error increases as the sampling 
interval changes from 2m to 5m. The reason for the error 
increase is that the search space for beacon node placement is 
pruned when increasing the sampling interval, or decreasing 
the selected area size. As discussed in Section Uvu we 
can take user-specified acceptable calculation time T acc and 
localization error A E acc as input to select a proper sampling 
interval and area. I.e., specifying T acc = 1632.23s, and 
e opt + A E acc < 2.85, we can select an area with size being 
^5 and interval being 2m for sampling. 


D. Comparison in Indoor Environment 

We apply four other placement methods, two random ones 
(Random and RKC) and two deterministic ones (Uniform and 
CERACC), as reference for comparison. The beacon nodes 
placed by these methods are shown in the sub-figures of Fig. 
[6] separately. The corresponding localization errors of these 
methods on EZPerfect and RADAR are shown in Table I VII and 
| VII| respectively. As can be seen from both of these two tables, 










































































































1) Random has the largest localization error; 2) the localization 
error CERACC is less than Random; 3) RKC and Uniform 
have roughly equal error, better than CERACC; 4) compared 
with the results of our placement method in Table [TV] and [VJ 
RKC and Uniform have approximately same performance in 
(^5, 5m) case, worse in other cases. 


E. Comparison in MetroFi 

We also compare placement methods in a dataset, MetroFi, 
of an outdoor environment. According to the previous result, 
we set T acc to 1631.23s (Table [IT} and 1229.61s (Table IP 
for EZPerfect and RADAR respectively, and let A E acc be 5m 
for both localization algorithms. Then, our method generates 
beacon node placement with the limitation of T acc and A E acc . 
The comparison results of EZPerfect and RADAR are shown 
in Table VIII| and [lx] respectively. As can be seen, 1) our 
placement method (BNP) has the lowest localization error; 2) 
Uniform is slightly better than RKC; 3) both Uniform and 
RKC have lower localization error compared with CEACC; 
4) Random has the worst performance. 


VI. Related Work 

Beacon node placement has been previously studied for 
coverage and distribution in research literature. In the coverage 
problem, it is required to use the minimal number of beacon 
nodes to achieving k-coverage in a given area. As it is 
proved to be NP-hard |9), some approximate (TJ and random 
CD placement methods are proposed to k-cover bounded or 
unbounded area. As for distribution, a typical work is to locate 
a target node by Trilateration or Multilateration when beacon 
nodes are in GDoP optimal distribution j26). And also, beacon 
nodes are deployed by random, max and grid placement 101 
for localization. 

Besides, localization algorithms also have inexplicit re¬ 
quirement for the coverage and distribution of beacon nodes. 
A considerable part of localization algorithms, such as 
Fingerprinting-based mmm and Proximity-based mm, 
require that beacon nodes should be placed to at least achieve 
1-coverage in the interest area, and prefer evenly scattered 
node placement . Some other algorithms such as Trilateration- 
based CUED require 3 or more coverage for solvability, and 
consider the distribution of beacon nodes with minimal GDoP 
integral value. Also, there exists algorithms such as MDS 
I23ll24l . SDP (6) and Hop-based |14][16|[Hl can work with 
few or even no beacon node (0-coverage), and the optimal 
beacon node distribution of these algorithms mainly depend 
on localization environment. 


VII. Conclusion and Future Work 

In this paper, we study the Beacon Node Placement (BNP) 
problem that beacon nodes should be deployed to minimize the 
localization error. We prove that BNP is NP-hard. In view of 
the hardness of BNP, we propose to synthesize a function /'(.) 
instead of /(.) to approximate the calculation of localization 
error. The main idea behind the synthesization process is 
to divide-and-conquer the calculation of /(.). We prove that 
the synthesization process executes with time complexity of 
0(T acc ) and generates the distribution of beacon nodes with 
localization error e — e opt < 7 A E acc . In the experiment, 
we test beacon node placement according to the generated 


TABLE VIII 

The localization error ol EZPerlect on MetroFi dataset 


EZPerfect 

ari. 

geo. 

med. 

abn. 

Random 

39.15 

32.73 

33.37 

26.38% 

RKC 

24.64 

19.83 

20.25 

4.16% 

Uniform 

23.02 

18.62 

19.83 

4.16% 

CERACC 

31.27 

26.27 

23.36 

11.11% 

BNP 

18.98 

16.67 

15.89 

1.38% 


TABLE IX 

The localization error of RADAR on MetroFi dataset 


RADAR 

ari. 

geo. 

med. 

abn. 

Random 

36.85 

29.46 

32.37 

23.61% 

RKC 

21.82 

17.27 

19.67 

2.77% 

Uniform 

20.52 

16.63 

18.26 

4.16% 

CERACC 

26.64 

23.83 

22.36 

9.72% 

BNP 

17.49 

15.74 

17.15 

1.38% 


distribution, and compare with other placement methods under 
various settings, such as 2910m 2 indoor floor, outdoor MetroFi 
dataset. The experimental results show the feasibility and 
effectiveness of our placement method. 

In the future, we will extend our placement method to han¬ 
dle the diversity of localization error, i.e., different localization 
accuracy may be wanted in different sub-area. This object can 
be achieved by introducing the weight model in localization 
error calculation. 
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